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ABSTRACT 

Dark matter or modifications of the Newtonian inverse-square law in the solar system are 
studied with accurate planetary astrometric data. From extra-perihelion precession and possi- 
ble changes in the third Kepler's law, we get an upper limit on the local dark matter density, 
Pdm ^ 3 x 10 -16 kg/m at the 2-a confidence level. Variations in the 1 jr 2 behavior are con- 
sidered in the form of either a possible Yukawa-like interaction or a modification of gravity 
of MOND type. Up to scales of 10 11 m, scale-dependent deviations in the gravitational accel- 
eration are really small. We examined the MOND interpolating function fi in the regime of 
strong gravity. Gradually varying fi suggested by fits of rotation curves are excluded, whereas 
the standard form fj,(x) — x/(l + x 2 ) 1 / 2 is still compatible with data. In combination with 
constraints from galactic rotation curves and theoretical considerations on the external field 
effect, the absence of any significant deviation from inverse square attraction in the solar sys- 
tem makes the range of acceptable interpolating functions significantly narrow. Future radio 
ranging observations of outer planets with an accuracy of few tenths of a meter could either 
give positive evidence of dark matter or disprove modifications of gravity. 
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1 INTRODUCTION 

Gravitational inverse-square law and its relativistic generalization 
have passed significant tests on very different length- and time- 
scales. Precision tests from laboratory and from measurements in 
the solar system and binary pulsars provide a quite impressive body 
of evidence, considering the extrapolation from the empirical basis 
I Adelb erser etaflbOOa IWillll200rt) . First incongruences seem to 
show up only on galactic scales with the observed discrepancy be- 
tween the Newtonian dynamical mass and the directly observable 
luminous mass and they are still in order for even larger gravita- 
tional systems. Two obvious explanations have been proposed: ei- 
ther large quantities of u nseen 'dark' m atter (DM) dominate the dy- 
namics of large systems fcwickvll933T) or gravity is not described 
by Newtonian theory on every scale iFinzil 19631) . Dark matter is all 
the general theory of relativity needs to overcome apparent short- 
comings and provides a coherent picture for gravitational phenom- 
ena from the laboratory to the cosmological context. The paradigm 
of cold DM when complemented with a positive cosmological con- 
stant (the so called ACDM scenario) is successful in explaining the 
whole range of galactic and extra-galactic body of evidence, from 
flat rotation curves in spiral galaxies to large scale structure forma- 
tion and evolution IPeacockl 19991) . 

The ACDM paradigm could be regarded as the definitive pic- 
ture apart from that the presumed existence of DM relies so long 

* E-mail: sereno@physik.unizh.ch (MS); jetzer@physik.unizh.ch (PhJ) 



only on its putative global gravitational effect, whereas direct de- 
tection by any independent mean is still lacking. This makes room 
to alternative proposals based on modifications of Newtonian grav- 
ity. In general, such proposals do not extend the inverse-square law 
to a regime in which it has never before been tested and they do 
not introduce any exotic component. Proposals are very different 
from each other. Some of them can make gravity stronger on scales 
of galaxies and explain flat rotation curves without dark matter 
(Milarom 1983); others realize a mechanism for the cosmic accel- 
eration without dark energy, for example as a result of gravity leak- 
ing on scales comparable to the horizon jDvali et aDl200Ct) . Two 
main alternative proposals have been discussed. In the first one, the 
gravitational potential deviates from the usual form at large dis- 
tances. A classical example is the inclusion of a Yukawa-like term 
in the gravitational potential. This is strictly related to more fun- 
damental theories where these additional contributions appear as 
the static limit of interactions due to the exchange of virtual mas- 
sive bosons l Adelberser et al. 2003). According to the second main 
choice, Newton's law fails when the gravitational acceleration is 
small rather than when the distance is large. The prototype and still 
one of the most empirically successful alternative to DM is Mil- 
grom's modified Newtonian dynamics (MOND) iMilgrom 1983; 
Sanders & McGauah 2002). With some basis in sensible physics, 
MOND can provide an efficient description of the phenomenol- 
ogy on scales ranging from dwarf spheroidal galaxies to cluster 
of galaxies b ut its cosmological extension is still in the childhood 
(Bek ensteinl2004) . 



2 M. Sereno & Ph. Jetzer 



High precision solar system tests could provide model inde- 
pendent constraints on possible modifications of Newtonian grav- 
ity. The solar system is the larger one with very well known mass 
distribution and can offer tight confirmations of Newtonian grav- 
ity and general relativity. Any deviation emerging from classical 
tests would give unique information either on dark matter and 
its supposed existence or on the nature of the deviation from the 
inver se-square law. Se veral authors have discussed this possibility. 
iTalmadge et alJ {l988) derived limits from the analysis of various 
planetary astrometric data set on the variation in the 1 /r 2 behav- 
ior of gravity. Experimental bounds on non luminous matter in so- 
lar orbit were derived either by considering the third Kepler's law 
lAnderson et alll989lll995h or by stu dying its effect upon perihe- 
lion precession ( Gron & Soleng 1996). The influence of a tidal field 
due to Galactic dark matter on the motion of the plane ts and satel- 
lites in the solar system was further i nvestigated bv lBraginskv et all 
119921) and lKlioner&Soffell ll993l) . The orbital motion of solar- 
system planets has been determined with higher and higher accu- 
racy lPitieval20055) and recent data allow to put interesting limits 
on very subtle e ffects, such as that of a non n ull cosmological con- 
stant Ijetzer & Serenoll2006tlSereno & JetzeJ2006t) . In this paper, 
we discuss what state-of-art ephemerides tell us about non Newto- 
nian or DM features. In section|2| we review standard expectations 
about Galactic dark matter at the solar circle and discuss some stan- 
dard frameworks for deviations from the inverse square-law, i.e. a 
Yukawa-like fifth force and the MOND formalism. Observational 
constraints from perihelion precessions and changes in the third 
Keplerian law are discussed in section [3] and [4] respectively. Sec- 
tionals devoted to some final considerations. 



2 BASICS 

Let us now briefly consider the main features of dark matter in the 
solar system and of some alternatives to Newtonian gravity. 



2.1 Dark matter 

In the dark matter scenario, Milky Way is supposed to be embedded 
in a massive dark halo. Realistic models of the Milky Way based 
on adiabatic compression of cold DM haloes can be b uilt in agree- 
ment with a full range of obs ervational constraints iKlypin et alJ 
l2002UCardone & Serenol2005l) . The local DM density at the solar 
circle is then expected to be pdm ~ 0.2 x 1CP 21 kg/m 3 , in ex- 
cess of nearly five orders of magnitude with respect to the mean 
cosmological dark matter density. 



2.2 MOND 

MOND underpins the principle that gravitation departs from New- 
tonian theory if dynamical accelerations are small. It was initially 
propo sed either as a modification of inertia or of gravity fMilgroml 
Il983l) . According to this second approach, the gravitational accel- 
eration g if related to the Newtonian gravitational acceleration gN 
as 



M(|g|/ao)g = gN 



(1) 



where ao is a physical parameter with units of acceleration and 
n(x) is an unspecified function which runs from fj,(x) = xatx<^i 
1 to fi(x) = 1 at x 3> 1. Whereas the Newtonian trend is recovered 
at large accelerations, in the low acceleration regime the effective 



gravitational acceleration becomes g ~ y/g^a . The asymptoti- 
cally flat rotation curves of spiral galaxies and the Tully-Fisher law 
are explained by such a modification and a wide range of obser- 
vations is fitted with the same value of ao — 1.2 x 10 -10 m s -2 
ISanders & McGaugj2002l) . 

The fx function is formally free but, as a matter of facts, fits to 
rotation curves or considera tions on the ext ernal field effects sug- 
gest a fairly narrow range l Zhao & Famaev 2006). The standard 
interpolating function proposed by Milgrom 1 1983), 



(j,(x) = xj\f\ 



+ x 2 



(2) 



provides a reasonable fit to rotation curves of a wide range of galax- 
ies. Based on a detailed stud y of the velocity curves o f the Milky 
Way and galaxy NGC 3198. iFamaev & Binnevj l2005h found out 
that interpolating functions which trigger a slower transition from 
the MONDian to the Newtonian regime should be preferred. They 
proposed the alternative interpolating function, 



fl(x) = x/(l + x). 



(3) 



Transition between the asymptotic regimes is smoother in Eq. J5} 
than in Eq. In principle, /i could be precisely determined from 
the observations of an ideal galaxy in which both the flat rotation 
curve and the luminosity distribution are known with high accuracy. 
The n function that best reproduces the Milky Way's rotation curve 
rotation seems to go smoothly fr om Eq. J2j at x < 1 to Eq. J3j at 
x > 10 jFamaev & Binnevl2 005). 

In the Newtonian regime, departures strongly depend on the 
way [i approaches 1 asymptotically. For a quite ge neral class of 
interpolating functions, we can write iMilgromligsl) 

yL{x) ~ 1 - k (l/x) m , (4) 

which leads to the modified gravitational field l Talmad ge et alJ 
Il988t) 

g~g N [l + fc (ao/|g N |) m ]. (5) 

For x ^> 1, Eq. J2j and Eq. J3} can be recovered for {ko,m} — 
{1/2, 2} and {1, 1}, respectively. 

Any viable relativistic theory embody ing th e MO ND 
paradigm, such as Bekentein 's TeVeS I Bek ensteinl 120041) or 
Sanders' BSTV l Sanders 2005), seems to require scalar and vecto- 
rial fields in addition to the usual tensor field. MOND phenomenol- 
ogy emerges as an effective fifth force associated with a scalar field. 
The interpolating function fj, is related to an auxiliary function of 
the scalar field strength. The parameterized post-Newtonian (PPN) 
formalism has been very effective in confronting metri c theories o f 
gravity with the results of solar-system experiments (Will 2006). 
Unfortunately, for the relativistic generalizations of MOND, the 
presence of both a scalar and a vector field, together with the free 
function in the Lagrangian that yields the expected dynamics in 
the low-acceleration limit, makes it problematic to derive the cor- 
responding PPN parameters. To date, preliminary derivations only 
concern the very inner solar system, where /i is very close to unity 
lBekensteirl2004l) . so that the very accurate determination of PPN 
parameters can not be directly used to test MOND. 



2.3 Yukawa-like fifth force 

Many long-range deviations can be characterized by an amplitude 
and a length scale. Let us consider additional contributions to the 
gravitational potential in the form of a Yukawa-like term, whose 
astrophysical consequences have been explored from the scale of 
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the solar system (Adelbergeretal. 2003) to the large scale structure 
of the universe jWhite & Kochanelj200AlAmendola & Ouercellini 
' 2004 ISealfon et al J 1200.4 IShirata et afll2005t ISereno & PeacocM 
2006). The weak field limit of the gravitational potential, cj), can 
be written as a sum of a Newtonian and a Yukawa-like potential; 
for a point mass M, 



1 + ay exp < — 



{-£} 



(6) 



where ay is a dimensionless strength parameter and Ay is a length 
cutoff. The potential in Eq. JSJ goes as oc 1/r both on a small 
scale (r <C Ay), with an effective coupling constant Goo (1 + ay), 
and on a very large scale, where the effective gravitational constant 
is Goo- We will take Goo = Gn/(1 + ay), so that the value of 
the coupling constant on a very small scale matches the observed 
laboratory value, Gn- The total gravitational acceleration felt by a 
planet embedded in the potential ^6) can be written as, 

:M 



~ Go 
-r — 



1 + ay 



(7) 



For ay < 0(> 0), gravity is enhanced (suppressed) on a large 
scale. The potential in Eq. |6| can be derived in a relativisti c grav- 
ity m odel that obeys the equivalence principle ( Zhvtni kov & Nesterl 
1994). A Yukawa-like contribution to the potential can be also con- 
nected to very-specific mass terms which appear in addition to the 
field theoretical analog of the usual Hilbert-Einstein Lagrangian 
( Babak & Grishchukl2003l) . 



3 PERIHELION PRECESSION 

As well known, a test body moving under the influence of the New- 
tonian potential of a central mass M will describe an ellipse with 
constant orbital elements. Due to a small, entirely radial perturba- 
tion, the arg ument of pericentre lu p will precess according to (see 
ISoffell 989. chapter 4) 



(1 



2Nl/2 



-8 An cos /, 



(8) 



where n = y GnM/o 3 is the mean motion of the unperturbed 
orbit, a the semimajor axis, e the eccentricity, / the true anomaly 
counted from the pericentre and SAr the radial component of the 
perturbing acceleration. The longitude of the ascending node is not 
affected. 

Data from space flights and modern astrometric methods made 
it possible to create very accurate planetary ephemerides and to 
precisely determine orbital elements of solar system planets. The 
latest EPM2004 ephemerides were based on more than 317000 po- 
sition observations collected over 1913-2003 and including radio- 
metric and optical astrometric observations of spacecraft, planets, 
and their satellites iPitieva 2005a). Such ephemerides were con- 
structed by simultaneous numerical integration of the equations of 
motion in the post-Newtonian approximation accounting for subtle 
effects such as the influence of 301 large asteroids and of the ring 
of small asteroids, as well as the solar oblateness. Extra-corrections 
to the known general relativistic predictions can be interpreted in 
terms of new physics. Results are listed in TableQ We considered 
the 2-cr upper bounds. When the additional non-Newtonian accel- 
eration is parameterized as constant, the average precession rate is 
given by 

=2U/2 



(1 



Best constraints on SAr come from Earth and Mars observations, 
sec TableQ 

Analyzed data from Pioneer spacecrafts cover an heliocentric 
distance out to ~ 70 AU and show an anomalous acceleration di- 
rected towards the Sun with a magnitude of ~ 9xl0~ 10 ms~ 2 
which first appeared at a distance of 20 AU from the Sun 
( Ande rson et alj Eo02). If such an acceleration were gravitational 
in origin it would be not universal. In fact, effects on orbits of in- 
ner and outer planets would be large enough to have been det ected 
given the present levels of accuracy ( And erson et aill2002l: lloricl 
2006b; ISandersl Eo06). The upper bound from Mars in TableQJis 
more than four orders of magnitude smaller then the Pioneer accel- 
eration. 



3.1 Dark matter 

Galactic dark matter can cause extra-perihelion precession in the 
solar system. A refined analysis should consider the anisotropy in 
the gravitational field in the solar system due to tidal forces induced 
by th e DM distribution JBraginskv et alJll992l : iKlioner & Soffell 
1993). As an alternative approach, a spheri cally symmetric distri- 
bution around the Sun can be considered iGron & Solendfl996t 
Khrip lovich & PitievalEoOr^ . In fact, the effect at a given orbital 
radius is essentially given by the total DM mass contained within 
it, w ith a very weak depe ndence on the actual density distribu- 
tion I Anderson et al. 1989). Dark matter density varies very slowly 
within the solar system and can be considered as nearly constant. 
Assuming a constant density pdm, the perturbing radial accelera- 
tion at a radius r is 8Ar = — (4-7rGNPDM/3)r. After substituting 
in Eq. {5J and averaging over a period, the extra-precession rate can 
be written as 



, . . 2Gn7TPDM /., 2\ V 2 

H) = U-e ■ 



(10) 



Note that for an effective uniform density of matter represented by 
a cosmological constant, i.e. pdm = — c 2 A/(47t Gn), the classical 
result for orbital prece ssion due to A is retrieved iKerr et alJ2003l : 
lletzer & Serenol2006l) . The best upper bound on local dark matter 
density comes from Mars data, see TableQ The accuracy on Mars 
precession should improve by more than six orders of magnitude to 
get constraints competitive with local estimates based on Galactic 
observables. 



3.2 MOND 

The rate of perihelion shift in the Newtonian regime of MOND 
(x 3> 1) with a generic interpolating function in the form of Eq. J4j 
can be expressed in terms of hypergeometric functions. Here, we 
report only the case of a small eccentricity when 



(wp) = —kon I — 



{1 + e 2 [l - m(5 - 2m)]/4 + C(e 4 )} , 



(11) 



-5 An. 



(9) 



where ryi = W Gs^M/ai). As for the DM case, the Mars data is 
the more effective in constraining the parameter space, see Fig.Q 
For ko ~ 1, we get m > 1.5. Results from solar system are 
in disagreement with expectations based on the extrapolation to 
the strong acceleration regime of the free functions preferred on a 
galactic dynamics basis. From the study of rotation curves, acceler- 
ations seems to continue to increase quite smoothly as 1 jr even in 
the intermediate MONDian regime. As a consequence, for x < 1 
free functions which trigger a smooth transition, as Eq. ^3), must 
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Table 1. 2-cr constraints from extra-precession of the inner planets of the solar system. 5w p is the observed extra-precession rate from Pitjeva (2005b); SAr 
is a constant perturbative radial acceleration at the planet orbit and pdm is the DM density within the planet orbit. 



Name 


Suj p (arcsec/year) 




Mb 


(m/ 


> 2 ) 






Pdm (kg/m 3 ) 


Mercury 


-0.36(50) X 10~ 4 


-1x10- 


12 < 


SA R 


< 


5x10- 


13 


< 4xl0- 14 


Venus 


0.53(30) x 10" 2 


-4x10" 


12 < 


SAr 


< 


6x10- 


11 


< SxlO" 14 


Earth 


-0.2(4) x 10~ 5 


-5x10- 


14 < 


SAr 


< 


3x10- 


14 


< 7xl0~ 16 


Mars 


0.1(5) x 10~ 5 


-3x10" 


14 < 


SAr 


< 


4x10- 


14 


< 3xl0" ie 



DO 
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Figure 1. Constraints on the MOND interpolating function, parameterized 
as in Eq. |4}, arising from extra-perihelion precession of inner planets. The 
shadow region is ruled out at the 2-cr confidence level. The points labeled 
Mil and F&B locate the interpolating function in Eq. |2} and 0, respec- 
tively. 



be preferred over expression s such as Eq. J2j l Famaev & Binned 
boost I Zhao & Famaev 2006). On the other hand, in the Newtonian 
regime, i> 1, the functional form of Eq. J3} is clearly excluded by 
solar system data, whereas Eq. J2j is still compatible. Combining 
data from solar system and galactic dynamics, in the comparison 
between Eq. J2j and Eq. 0, the first one seems to be preferred 
in both the deep MONDian and Newtonian regimes, whereas the 
second one gives a be tter fit for the i ntermediate region. Similar 
considerations induced Is andersl |2006) to argue that the total grav- 
itational acceleration is strictly Newtonian, i.e. oc 1/r 2 , on small 
scales and that the transition to the total asymptotic acceleration 
oc 1/r shows up through a plateau region between 10 2 and 10 3 AU 
where the extra acceleration is more or less constant. The accuracy 
on Mars data should improve by nearly four orders of magnitude to 
disprove the standard interpolating function in Eq. J2j. 

3.3 Yukawa fifth force 

The anomalous precession rate due to a Yukawa-like contribution 
to the gravitational potential is 



(12) 



Figure 2. Constraints on the Yukawa-like fifth force parameters, for positive 
«Y, arising from extra-perihelion precession of inner planets. The shadow 
region is ruled out at the 2-cr confidence level. 
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Extra-precession data for a planet with semimajor axis a mainly 
probe scale lengths of Ay ~ a/2. Solar system data allow to 
constrain departures from the inverse-square law with high accu- 
racy fora scale length Ay ~ 10 10 — 10 11 m l Talmad ge et alJ 



1988; Iorio 2005). Bounds are mainly determined from Mercury 



and Earth data, see Figs. [2] and [3] For Ay 

-5 x 10 -11 <q y <6x 10 -11 . 



10 m, we get 



4 THIRD KEPLER'S LAW 

A departure from the inverse-square law could affect the radial mo- 
tion of a body around a central mass and a change in the Kepler's 
third law would occur. The Newtonian law of motion for a test body 
in a circular orbit around a central mass M, in presence of a per- 
turbing radial acceleration SAr, can be written as 



2 

lu r 



GnM 



5A T 



(13) 
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Figure 3. Constraints on the Yukawa-like fifth force parameters, for neg- 
ative c*y. arising from extra-perihelion precession of inner planets. The 
shadow region is ruled out at the 2-cr confidence level. 



Table 2. 2-cr upper bounds from anomalous mean motion of the solar 
system planets; 8a is the uncertainty on the semimajor axis from Pitjeva 
(2005a); 8Au is an anomalous constant radial acceleration; pdm is the 
dark matter density. 



competitive with data from perihelion precession. Errors in Table|2| 
are formal and could be underestimated. Current accuracy can be 
determined evaluating the discrepancies in different ephemerides. 
Differences in the heliocentric distances do not exceed 10 km for 
Jupiter and amount to 180, 410, 1200 an d 14000 km fo r Saturn, 
Uranus, Neptune and Pluto, respectively lPitievall2005bl) . Bounds 
from outer planets reported in Table [2] should be accordingly in- 
creased. 

4.1 Dark matter 

Bounds on pom from deviations in the mean motion of inner 
planets, see Table [2] are of the same order of magnitude of con- 
straints from extra-precession. Observations of outer planets pro- 
vide constraints that are an order of magnitude larger but they 
give the best future prospects. Unlike inner planets, radio-technical 
observations of outer planet are still missing and their orbits can 
not be determined with great accuracy. Since the required accu- 
racy to probe the effects of a given uniform background decreases 
as oc a~ 4 , whereas the measurements precision of ranging ob- 
servations is roughly proportional to the range distance, explo- 
ration of outer planets seems pretty interesting. Dark matter with 
Pdm — 0.2 x 10~ 21 kg/m 3 could be detected if the orbital axis of 
the Uranus, Neptune and Pluto orbits were determined with an ac- 
curacy of 5a ~ 3 x 10~ 2 , 2 x 10" 1 and 5 x 10" 1 m, respectively. 
Up till now, the only ranging measurements available for Uranus 
and Neptune are the Voyager 2 flyby data, with an accurac y in the 
determination of distance of ~ 1 km (And erson et alll9 95). not so 
far from what required to probe solar system effects of dark mat- 
ter. 1 



Name 


8a (m) 


\SAn\ O/ 


; 2 ) 


PDM (kg/m 


Mercury 


0.105 x 10+° 


< 4xl0~ 


13 


< 3x10- 


14 


Venus 


0.329 x 10+° 


< 2x10" 


13 


< 7x10- 


15 


Earth 


0.146 x 10+° 


< 3xl0~ 


14 


< 8x10- 


16 


Mars 


0.657 x 10+° 


< 4x10" 


14 


< 7x10- 


l(i 


Jupiter 


0.639 x 10+ 3 


< lxlO" 


12 


< 5x10- 


16 


Saturn 


0.4222 x 10+ 4 


< lxlO" 


12 


< 3x10- 


15 


Uranus 


0.38484 x 10+ 5 


< 1x10- 


12 


< 2x10- 


15 


Neptune 


0.478532 x 10+ 6 


< 4x10- 


12 


< 3x10- 


15 


Pluto 


0.3463309 x 10+ 7 


< 1x10- 


11 


< 8x10- 


15 



4.2 MOND 

Results from analysis of mean motion are similar to extra- 
precession analysis. The interpolating function in Eq. |3) is not con- 
sistent with solar system data. From Uranus data, we get m > 1.4 
assuming ko — 1. Again, best future prospects are related to radio- 
technical determination of orbits of outer planets. The standard in- 
terpolating function in Eq. could be (dis-)probed if the axes of 
the Uranus, Neptune and Pluto orbit were determined with an ac- 
curacy of 5a ~ 3 x 10, 3 x 10 2 and 1 x 10 3 m, respectively. 



where uj is the angular frequency and M e s = M(l + <5_4 r /.4N) 
is the effective mass felt by the orbiting planet. In other words, 
the angular frequency will differ from the mean motion n = 
^GnM/a 3 . It is 

Sn _ 1 5A T 
n ~ 2 .4 N ' 



(14) 



Variation of the effective solar mass felt by the solar system in- 
ner planets with respect to the effective ma sses felt by oute r plan- 
ets could probe new physics lAnderson et alj|l989L Il995l) . We 
can evaluate the statistical error on the mean motion for each 
major planet from the uncertainty on the semimajor axis, 5n — 
— (3/2)nSa/a, and translate it into an uncertainty on the effec- 
tive acceleration. Results for a constant additional acceleration term 
are listed in Table [2] Assuming 5A r being constant, it would be 
< 5 x 10" 12 m/s 2 in the range 20-30 AU, as can be inferred 
from Uranus and Neptune orbi ts, well below the anomalous Pio- 
neer acceleration llorio 2006b). Limits from Earth and Mars are 



4.3 Yukawa fifth force 

Comparison of Keplerian mean motions of inner and outer planets 
can probe a Yukawa-like contribution only if planets feel different 
effective gravitational constants. Such test is insensitive to values 
of Ay either much less the orbit radius of the inner planet o r muc h 
larger than the orbit of the outer planets l Adel berger et alj|2003l) . 
Differently from extra-precession of perihelion, which appears only 
for departures from the inverse square law, changes in the mean 
motion can appear even if both planets feel a gravitational accel- 
eration oc 1 /r 2 but with different renormalized gravitational con- 
stants. Considering inner planets, Earth data give | ay | < 6 x 10 _ 12 



1 The approach followed here is more conservative than a similar analysis 
appeared in Iloriol f2()06a). where it is assumed that all of the observational 
residuals can be fully accounted for with a suitable combination of different 
effects. Furthermore, lloriol {2006a ) considers a quite peculiar form of the 
gravitational potential valid only for a dark mass distribution fully contained 
within the considered orbit. 
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Figure 4. Constraints on Yukawa fifth force parameters, for the absolute 
value of »y, arising from deviations from the third Kepler's law. The 
shadow region is ruled out at the 2-cr confidence level. 

for Ay 5 2 x 10 10 m. The best constraint from outer planets is due 
to Jupiter, with \a Y \ < 5 x 10~ 9 for A Y < 10 11 m. 



5 CONCLUSIONS 

Debate between dark matter and departures from inverse-square 
law is still open. Considering both theoretical and observational 
aspects, dark matter seems to be slightly preferred. If on a galac- 
tic scale the two hypotheses match, on the cosmological side only 
DM can give a consistent framework. This might shortly change 
with the steady improvements in relativistic generalization of the 
MONDian paradigm. So, in our opinion, it is of interest to exam- 
ine results on a very different scale, that of the solar system. Solar 
system data have been confirming predictions from the general the- 
ory of relativity without any need for dark matter and it is usually 
assumed that deviations can show up only on a larger scale. In this 
paper, we have explored what we can learn from orbital motion of 
major planets in the solar system. Results are still non-conclusive 
but nevertheless interesting. Best constraints come from perihelion 
precession of Earth and Mars, with similar results from modifica- 
tions of the third Kepler's law. The upper bound on the local dark 
matter density, pdm $ 3xl0 -16 kg/m 3 , falls short to estimates 
from Galactic dynamics by six orders of magnitude. 

Deviations of the gravitational acceleration from 1/r 2 are re- 
ally negligible in the inner regions. A Yukawa-like fifth force is 
strongly constrained on the scale of ~ 1 AU. For a scale-length 
Ay ~ 10 11 m, a Yukawa-like modification can contribute to the to- 
tal gravitational action for less then one part on 10 11 . Similar limits 
could be achieved by precise measurements on the proof masses 
carried on board of the LISA Pathfinder satellite 2 (Speake, private 
communication). In fact, instantaneous measurements of the drag- 
free test-mass acceleration during the transfer orbit towards the first 

2 http://www.rssd.esa.int/index.php?project=LISAPATHFINDER 



Sun-Earth Lagrange point could in principle test the inverse square 
law on a scale length of ~ 1 AU (Speake, private communication). 
Results on a similar scale-length could be obtained through a de- 
tailed analysis of binary pulsars. The periastron shift, the gravita- 
tional redsfhift/second-order Doppler shift parameter and the rate 
of change of orbital period are sensitive to scalar-tensor gravit y and 
to any other deviation from the general theory of relativity IWilll 
2006). Dipole gravitational radiation associated with violations of 
the equivalence principle in its strong version could cause an ad- 
ditional form of gravitational damping and a significant change of 
the orbital period could occur, in particular for a binary pulsar sys- 
tem with objects of very dissimilar mass l Will 2006). A massive 
graviton associated with a Yukawa-like fifth force could also affect 
the speed of propagation of gravitational waves and induce r adia- 
tion e ffects at the reach of future gravitational wave detectors IWilll 
2006). 

A large class of MOND interpolating function is excluded by 
data in the regime of strong gravity. The onset of the asymptotic 
1 jr acceleration should occur quite sharply at the edge of the so- 
lar system, excluding the more gradually varying \i(x) suggested 
by fits of rotation curves. On the other hand, the standard MOND 
interpolating function (i(x) = x/(l + x 2 ) 1 ^ 2 is still in place. Stud- 
ies on planetary orbits could be complemented with independent 
observations in the solar system. Mild or even strong MOND be- 
havior might become evident near saddle points of the total gravi- 
tational potential, where MONDian phenomena might be put at the 
reach of mea surements by spacecraft equipped with sensitive ac- 
celerometers iBekenst ein & Masueiio 2006). As a matter of fact, 
fits to galactic rotation curves, theoretical considerations on the ex- 
ternal field effects and solar system data could determine the shape 
of the interpolating function with a good accuracy on a pretty large 
intermediate range between the deep Newtonian and MONDian 
asymptotic behaviors. 

Future experiments performing radio ranging observations of 
outer planets could greatly improve our knowledge about gravity 
in the regime of large accelerations. The presence of dark matter 
could be detected with a viable accuracy of few tenths of a meter 
on the measurements of the orbits of Neptune or Pluto, whereas 
an uncertainty as large as hundreds of meters would be enough to 
disprove some pretty popular MOND interpolating functions. 

In order to become really competitive with general relativ- 
ity and the ACDM paradigm, MOND should be predictive on the 
whole range of observed systems from solar system to the cos- 
mic microwave background radiation. On a galactic scale, effects 
of DM or MOND are pretty similar and very difficult to distin- 
guish each other but there might be some detectable differences 
on a smaller scale. In fact, the local value of DM at the solar cir- 
cle is pretty much fixed by Galactic dynamics whereas the MOND 
behavior in the regime of strong accelerations probed locally is 
not univocal on a theoretical and observational basis. Neverthe- 
less, only a very small class of interpolating free functions would 
give the same perturbation on the orbits of outer planets as that 
from local DM. Matching the expectations from DM with future 
radio ranging observations would be an important, nearly conclu- 
sive confirmation of its existence. On the other hand, deviations at 
a different order of magnitude, as those expected for a large variety 
of MOND interpolating functions, would be a strong indication of 
departure from the inverse-square gravitational law. 
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